Superscaling in Dilute Fermi Gas and Its Relation to General Properties of 
the Nucleon Momentum Distribution in Nuclei 



O 

o 

(N 



A.N. Antonov, 1 M.V. Ivanov, 1 M.K. Gaidarov, 1 ' 2 and E. Moya de Guerra 2 - 3 

'institute for Nuclear Research and Nuclear Energy, 
Bulgarian Academy of Sciences, Sofia 1784, Bulgaria 
2 Instituto de Estructura de la Materia, CSIC, Serrano 123, E-28006 Madrid, Spain 
3 Departamento de Fisica Atomica, Molecular y Nuclear, Facultad de Ciencias Fisicas, 
Universidad Complutense de Madrid, E-28040 Madrid, Spain 

The superscaling observed in inclusive electron scattering is described within the dilute Fermi gas 
model with interaction between the particles. The comparison with the relativistic Fermi gas (RFG) 
model without interaction shows an improvement in the explanation of the scaling function f(ip') 
in the region ip' < — 1, where the RFG result is f(ip') = 0. It is found that the behavior of f(ip') 
for ip' < — 1 depends on the particular form of the general power-law asymptotics of the momentum 
distribution n(k) ~ l/fc 4+m at large k. The best agreement with the empirical scaling function is 
found for m ~ 4.5 in agreement with the asymptotics of n(k) in the coherent density fluctuation 
model where m — 4. Thus, superscaling gives information about the asymptotics of n(k) and the 
NN forces. 

PACS numbers: 25.30.Fj, 21.60.-n, 21.10.Ft, 21.10.Gv 
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The concepts of y-scaling 0, i, i, i, i, i, 0, [EJ 
and superscaling (based on ^'-scaling variable) [lC 
El EE EI EE EE EE El] have been used for extensive 
analyses of the vast amount of inclusive electron scatter- 
ing world data (see also [l9|]). These analyses showed 
the existence of high-momentum components in the nu- 
cleon momentum distribution (MD) n(k) at momenta 
k > 2 fm _1 due to the presence of nucleon- nucleon (NN) 
correlations. Scaling of the first kind (i.e. no dependence 
on the momentum transfer) can be observed at excita- 
tion energies below the quasielastic (QE) peak. Scaling 
of second kind (i.e. no dependence on the mass number) 
is excellent in the same region. When scaling of first and 
second type occur one sots that superscaling takes place. 
It was shown (e.g., in [IE EE E3 Ell) that the physi- 
cal reason for the superscaling phenomena is the specific 
high-momentum tail of n(k) which is similar for all nu- 
clei. 

It has been shown [IE HH that due to the contribution 
introduced by inelastic scattering together with the corre- 
lation contribution and meson exchange currents [IE HH 
both scaling of the first and, to a lesser extent, of the 
second kind are violated at energies above the QE peak. 

The theoretical concept of superscaling has been intro- 
duced in [IE EH using the properties of the relativistic 
Fermi gas (RFG) model . The Fermi momentum for the 
RFG was used as a physical scale to define the proper 
scaling variable ip' for each nucleus. As emphasized 
in [IE], however, the actual dynamical physics content of 
the superscaling phenomenon is more complex than that 
provided by the RFG model. In particular, the extension 
of the superscaling property to large negative values of 
ip' (ip' < -1) is not predicted by the RFG model. The 
QE scaling function in the RFG model f^p G (ip') = 
for ip' < —1, whereas the experimental scaling function 



/Q E (?//) has been observed for large negative values of 
ip' up to ip' w —2 in the data for (e, e') processes. Thus, 
it has been necessary to consider the superscaling in the- 
oretical approaches which go beyond the RFG model. 
One of them is the coherent density fluctuation model 
(CDFM) (e.g., 0, HH) which gives a natural extension 
of the Fermi gas case to realistic finite nuclear systems. It 
was shown in [lj, EE EE El that in the CDFM both ba- 
sic quantities , density and momentum distributions, are 
responsible for the scaling and superscaling phenomena 
in nuclei. The results of the CDFM for the QE scaling 
function f(ip') agree with the available experimental data 
at different transferred momenta and energies below the 
QE peak position, showing superscaling for ip' < 0, in- 
cluding ip' < — 1 and going well beyond the RFG model. 
Secondly, as pointed out in [IE EH , the nucleon momen- 
tum distribution for various nuclei obtained in [26| (and 
with the modification in [Hj]) within a theoretical ap- 
proach based on the light-front dynamics (LFD) method 
(e.g., [13, HI) can also be used to describe both y- and 
^'-scaling data. 

The superscaling analyses of inclusive electron scatter- 
ing from nuclei (for energies of several hundred MeV to a 
few GeV) have been extended in 29] to include not only 
QE processes but also the region where A excitation dom- 
inates. Both QE scaling functions fQ E (ip') and / a (V'a) 
in the delta region were deduced from phenomenological 
fits to electron-nuclei scattering data. Generally, the spe- 
cific features of the scaling function should be accounted 
for by reliable microscopic calculations that take final- 
state interactions into account. In particular, the scal- 
ing function f® E (ip') with asymmetric shape obtained 
in [IE HH by using a relativistic mean field (RMF) for 
the final states agrees well with the experimental scaling 
function. 

The features of the superscaling phenomenon in inelas- 
tic electron scattering have induced studies of neutrino 
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scattering from nuclei on the same basis. The neutrino- 
nucleus interactions have been studied using the super- 
scaling analyses of few-GeV inclusive electron scatter- 
ing data in a method proposed in [29] to predict the 
inclusive vA and vA cross sections for the case of 12 C 
in the nuclear resonance region. Various other theoret- 
ical considerations {e.g., [II IE [H, HE [H, 13, H [H, 
liol . Ull . S3, |43|) h ave been devoted to studies of both 
neutral- (e.g. j3|j3l^j3E |3j) and charge-changing 
{e.g., [H, [H M, M, M, IS mL S3 S3, El Hi) neutrino- 
nucleus scattering. 



The CDFM and the LFD method have been extended 
in [l8[ from the QE to the delta-excitation region of the 
inclusive electron scattering and the QE scaling func- 
tions calculated in both methods have been used to cal- 
culate and to predict charge-changing neutrino-nucleus 
cross sections of the and (77 M ,^ + ) reactions on 

12 C at energies from 1 to 2 GeV. The asymmetry in the 
CDFM QE scaling function has been introduced in a phe- 
nomenological way. These analyses make it possible to 
gain information about the nucleon correlation effects on 
both nucleon momentum and local density distributions. 
It became clear that only the detailed knowledge of the 
behavior of n{k) at high momenta in realistic nuclear 
systems could lead to quantitative agreement with the 
experimental scaling function. On the other hand, the 
behavior of the latter gives a valuable information about 
the NN correlation effects on the tail of the momentum 
distribution. So, it was shown in [l6| within the CDFM 
that the y- and ^'-scaling data are informative for n{k) 
at momenta up to k ps 2-2.5 fm -1 and it was concluded 
that further experiments are necessary in studies of the 
high-momentum components of the nucleon momentum 
distributions. 



II. THEORETICAL SCHEME AND RESULTS 
OF CALCULATIONS 



In the first part of this Section we consider the hard- 
sphere dilute Fermi gas, i.e. the low-density Fermi gas 
whose particles interact via a repulsive hard-core poten- 
tial (see, e.g., 0, EE SI, S3, El] ) , and use the nucleon 
momentum distribution in such a system to calculate 
the scaling function / HSDFG (-(//). The quantities of in- 
terest in the HSDFG model can be expanded in powers 
of the parameter kpc, where c (> 0) denotes the hard- 
core radius in the NN interactions or it is identified with 
the scattering length in free space, and kp is the Fermi 
momentum. In [44j the value of fc^c is adopted to be 
equal to 0.70 which corresponds to NN core radius of 
c = 0.50 fm and the typical value of the Fermi momentum 
kp = 1.40 fm -1 . As was pointed out by Migdal [iq] . n{k) 
in the normal Fermi gas is discontinuous at the Fermi 
momentum. The analytical expressions for the dimen- 
sionless n{k) in the HSDFG obtained in [H,E1| have the 
form: 



n{k) = n < {k)+n > {k) with 



7i<(fc) = for k>kf 
n> {k) — for k<kp 



(1) 



At fc < kF'- 



n < {k)^l-^-l{k F c) 2 



(71n2-8)x 3 + (10-31n2)x 



, 1+x , {2-x 2 ) 1 / 2 +x 

+2 In 2{2-x 2 f' 2 In ' 

1—X {2 — x z ) L ' z — x 

where x = k/kp and v — 4 [44| is adopted. 
At 1 < x < V2: 



(2) 



The aim of the present work is to consider in more 
detail the connection between the NN forces in nuclear 
medium and their effect on the components of n{k) from 
one side and, from the other side, the role of n{k) on the 
behavior of the QE scaling function. For this purpose 
we use firstly the MD in a hard-sphere dilute Fermi gas 
model (HSDFG) {e.g., @, EE EE S3, El) to calculate the 
scaling function. Secondly, we make an attempt to throw 
light on the connection between the generally established 
high-momentum asymptotics of n{k) 0, EE El, HI, HH 
and the QE scaling function. The latter makes it pos- 
sible to establish (at least approximately) the particular 
form of the power-law decrease of n{k) at large values 
of k. This makes it possible to extract additional infor- 
mation about the NN forces from the description of the 
superscaling phenomenon. 



The theoretical scheme and the results of calculations 
are given in Section ITT] The conclusions are summarized 
m Section HTT1 
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FIG. 1: The momentum distribution n(k) in hard-sphere di- 
lute Fermi gas [44J as a function of x — k/kF- 
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FIG. 2: The scaling function f(t[>') in HSDFG calculated for 
different values of &fc in comparison with the RFG model re- 
sult. The grey area shows experimental data taken from [l3| ]. 



At x > 3: 



n>(fc)=2-^-(fc F c)^21n 
Sir z x 
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The momentum distribution in the HSDFG model is 
presented in Fig. [1] for the value of kpc = 0.70. 

Following the definition of the ^'-scaling function given 
by Barbaro et al. [llj , one can obtain for the case of the 
HSDFG: 

OO 

/H BD ro( ^ )= 3 J (6) 

ICI/V 

where r/p = hp /win, itin being the nucleon mass, and 

C = i>' { [ yjl + if F l] [2 + r ( yjl + Vp l) ] } • (7) 

In Eqs. ((6]) and ([7|) the dimensionless scaling variable 
ip' 2 (in units of the Fermi energy) has the physical mean- 
ing of the smallest kinetic energy which one of the nuclc- 
ons responding to an external probe can have . Since 
rip <C 1 we write Eq. © as: 



i) the Fermi momentum distribution of the initial nu- 

3 

cleon in the nucleus ^-9(kp — \k\) is replaced 

by the MD (with dimension), namely J P(k, E)dE, 
where P(k, E) is the spectral function, and 

ii) for k > kp the step function 8(kp — \k\) is retained 
to take into account approximately Pauli blocking 
for the final nucleon. 

In Figs. [2] and [3] we give the results for the HSDFG scal- 
ing function (in logarithmic and linear scale, respectively) 
calculated for different values of kpc from 0.70 to 0.28 
and compared to the result for the scaling function in 
the RFG model. One can see that the HSDFG scaling 
function is extended for large negative values of ip' in con- 
trast to the case of the RFG scaling function, but there is 
not a good agreement with the experimental data. One 
can see in the figures also the step behavior of the scal- 
ing function which reflects the discontinuity of n(k) at 
k = kp. Due to these reasons we consider in the second 
part of this Section the relation between the asymptotic 
behavior of the MD and the -(//-scaling function. 

A relationship between the NN force in nuclear 
medium Vnn(^) and the asymptotic behavior of the mo- 
mentum distribution n(k) was derived with great general- 
ity in 49, 50, Hlj. In these articles Amado and Woloshyn 
showed that the asymptotics of n{k) at large values of k 
is a power-law decrease: 



t HSDFG 



xn(x)dx. 



(8) 



It can be seen from Eq. (|SJ) that, as expected, the HSDFG 
system under consideration also exhibits superscaling. 
Eq. ([SJ) has been obtained using the approximations fol- 
lowing [431 ] : 



n(k) 



V NN (k) 



k 2 



(9) 



where Vnn(&) is the Fourier transform of the NN inter- 
action VnnM- I n the case of delta- forces (i.e. in the 
HSDFG) the asymptotics is n(k) - 1/fc 4 45]. It is still 
unknown if k or k/A must be large for Eq. @ to ap- 
ply [49|, [5^, HH- I n principle, it is shown in |44| that 
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FIG. 3: The same as in Fig. [2] but a linear scale for /(■(/>') 




FIG. 4: The momentum distribution in HSDFG n(x) multi- 
plied by x 4 = (k/kp) 4 . 



n(k) in the HSDFG decreases faster than k 4 , typically 
like ~ l/fc 4+m , where m > 0. In Fig. H we show x n(x) 
as a function of x — k/kp. It can be seen that the HS- 
DFG n(k) decreases approximately like ~ l/k A+m with a 
small value of to, because the result for x n(x) is almost 
a straight line which decreases slowly with the increase 
of x = k/kp > 1. 

Further in this work we study the question about the 
general feature of the NN force Vnn(&) that results in an 
n(k) with a power-law behavior that best agrees with the 
scaling function. To this aim we assume an NN-potential 
Vnn(^) different from a delta- function and calculate the 
scaling function using different asymptotics for n(k) in 
the dilute Fermi gas at k > kp. Therefore we look for 
the proper value of to. For k < kp we use n(k) [Eq. (2)] 
from [441 . but for k > kp we use: 




FIG. 5: The scaling function in a dilute Fermi gas calculated 
using Eq. (|13|l for different values of m in the asymptotics of 
the momentum distribution n(k) ~ l/k 4+m given in compar- 
ison with the RFG result. The grey area shows experimental 
data taken from [l3| . 



of n{k) and is equal to 



^ = 5|!(l + m )4+ m . 



(11) 



The factor 0.24 corresponds to the result for the part of 
the normalization (for k > kp) from the total normaliza- 
tion condition: 



n{k)d 3 k = 1. 



(12) 



Finally, from Eq. ([8]) one can obtain the following ex- 
pression for the scaling function: 



/(V0=0.T2 



1 



1 



l\2+m ' 



(13) 



In Fig. [5] we present the results for the scaling function 
[Eq. ([5])] for different values of to, compared with the 
RFG model result. One can see that agreement with the 
experimental QE scaling function is achieved when the 
value to sa 4.5 is used in Eqs. UHl), d) and (p3|). This 
means that the power-law decrease of n(k) which gives 
an optimal agreement with the data is 



n{k) 



1 



(14) 



We should note that this particular form of the 
power-law asymptotics is close to that obtained in the 
CDFM [H 



n(k) 



1 



(15) 



n(k) = N 



1 



k 4+r, 



for k > kp. 



(10) 



i.e. it corresponds to n(k) ~ l/k 4+m with to = 4. The 
inverse Fourier transform of Vnn(^) for to = 4 and to = 5 
The value of N is obtained by the total normalization gives T^nn(^) ~ 1/r and Vnn(^) ~ 1/r 1 / 2 , respectively. 
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We would like to emphasize the consistency of both 
the optimal asymptotics of n(k) for the dilute Fermi gas 
found in this work [Eq. ([Til)] with that in the CDFM 
[Eq. (f~~)]. As was shown in~fla.fTl.IT7l.IT8j the calculated 
QE scaling function f(ip') in the CDFM agrees well with 
the experimental scaling function. This fact shows that 
the behavior of the QE scaling function depends mainly 
on the particular form of the power-law asymptotics of 
the nucleon momentum distribution. This is proved on 
our work by the similarities of the result for the case of 
an interacting dilute Fermi gas with that one obtained in 
the CDFM as a model accounting for NN correlations in 
realistic finite nuclear systems. 

III. CONCLUSIONS 

The results of the present work can be summarized as 
follows: 

i) The superscaling observed in inclusive electron scat- 

tering from nuclei is considered within the model 
of dilute Fermi gas with interactions between par- 
ticles. The latter gives an improvement in com- 
parison with the results of the relativistic non- 
interacting Fermi gas model allowing one to de- 
scribe the QE scaling function for t/j' < —1, whereas 
the RFG model gives f(i/j') = in this region. 

ii) It is established that the hard-sphere (with delta- 

forces between nucleons) approximation for the di- 
lute Fermi gas is quite a rough one. The use of 
more realistic NN forces leading to to ~ 4.5 in- 
stead of m — (for delta-force) in the well-known 
power-law asymptotics of the momentum distribu- 
tion n{k) ~ l/fc 4+m at large k leads to a good ex- 
planation of the data for the ^'-scaling function in 
inclusive electron scattering from a wide range of 
nuclei. 

iii) The asymptotics of n(k) ~ 1/fc 8 5 found in the dilute 

Fermi gas by optimal fit to the data f or f (ijj') is 
similar to that in the CDFM (~ 1/fc 8 ) [2jf which, 
being a theoretical correlation model, describes the 
superscaling in the quasielastic part of the electron- 
nucleus scattering. Thus, the momentum distribu- 
tion in the dilute Fermi-gas model with realistic NN 



forces can serve as an "effective" momentum distri- 
bution (a step-like one with a discontinuity) which 
gives a similar result for f(ip') as the correlation 
methods for realistic finite nuclear systems. It can 
be concluded that the momentum distribution with 
asymptotics from ~ 1/k 8 to ~ 1/k 8 ' 5 is the proper 
one to explain the phenomcnological shape of the 
scaling function obtained from inclusive QE elec- 
tron scattering. 

As already mentioned, the superscaling is due to the 
specific high-momentum tail of n(k) similar for all nuclei 
which is known to be caused by the short-range and ten- 
sor correlations related to peculiarities of the NN forces 
near their core. The main result of the present work 
might be the observation that the values of f(ip') for 
tp' < — 1 depend on the particular form of the power-law 
asymptotics of n(k) at large k which is related to a corre- 
sponding particular behavior of the in-medium NN forces 
around the core. Namely, we point out that the power- 
law decrease of n(k) as ~ l/fc 4+m with m ~ 4.5 in the 
interacting dilute Fermi gas is the proper one and it is 
close to that obtained in CDFM (m = 4 [25]) which de- 
scribes the superscaling correctly as well. The NN force 
for m — 4 is expected to go as Vnn(^) ~ 1/r and for 
rn = 5 to go as Vnn(0 ~ (1/r) 1 / 2 . Hence, the present 
study allows one to conclude that the important property 
of the repulsive short-range core (leading to NN correla- 
tions and high- momentum tail of n{k)) is that it goes to 
infinity for r — > as 1/r or softer. The link between 
the asymptotic behavior of n{k) and NN force implies 
that inclusive QE electron scattering from nuclei provides 
important information on the NN forces in the nuclear 
medium. 
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